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ABSTRACT 

A narrow band wave spectrum was applied to theoretical 
relationships previously developed for set down and set up 
in an attempt tofind second order non-steady solutions 
for these concepts. The initial effort was to apply this 
spectrum to the radiation stress tensor using linear wave 
theory. Another development was attempted by incorporating 
the spectrum and the solution of the long wave equation in- 
to the Bernoulli and vertical momentum equations. Results 
obtained indicate that the solution for mean water level 
outside the surf zone is composed of a steady component and 
a periodic unsteady component; the periodic component being 
of the form of a long wave with a frequency lower than the 
components of the wave spectrum. The solution for set up 
is then composed of the same type components. The exact 
relationships depend on the patching process that is made 
for the solutions through the breaker line. 
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I . I NTRODUCT I ON 



Most of the observable phenomena along a coastline are 
the direct result of the action of the incoming waves, 
waves which begin in most instances as a disorganized con- 
fused state of the ocean surface, produced by a storm far 
at sea. In their transit across the vast expanse of the 
ocean, they begin to sort themselves out and form a somewhat 
regular oscillation of the ocean surface, the lower frequen- 
cy oscillations traveling faster and thereby leading the 
train. This train eventually ends by encountering a beach, 
where its energy is expended in the form of breaking waves. 

Since it is this aspect of the wave's life cycle that 
influences man the most, considerable effort has been ex- 
pended investigating this area. Among other effects, it has 
been observed that waves, in the process of shoaling and 
eventual breaking, produce a variation in the mean sea level 
This variation in sea level has been considered as the pri- 
mary cause for such nearshore currents as rip currents. The 
sea surface variation consists of: (a) a gradual depression 

of the mean sea level beginning offshore and reaching a maxi 
mum at the breaker line and (b) inside the surf zone a slooe 
of mean sea level which increases and extends shoreward to 
a point on the beach higher than the still water line. The 
depression is termed set down and the slope is called set up 
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It is reasonable to expect that variations in the amount of 
set down or set up along a beach can provide the head to 
produce a current. 

Previous investigations have considered steady state 
solutions. Experimental results obtained by both Bowen 
Cl967l] and Van Dorn Cl976l] agree quite favorably with the 
steady state solutions produced by Lo ng ue t-H i gg i n s and 
Stewart Cl962j using linear wave theory. 

This investigation considers an application of a simple 
wave spectrum to the existing theories in an attempt to ob- 
tain a non-steady solution for the set down and set up 
phenomena. An opening chapter on background is provided to 
ensure the necessary understanding of the existing theories 
and their development. Included is a section on the develop- 
ment of the "radiation stress tensor", a concept which was 
proved useful by Longuet-H i gg i ns and Stewart Cl962j in 
treating the shoaling process of waves. 

Chapter III deals with the application of a simple 
spectrum to the derivation of set down and set up. First 
order linear theory and first order theory including a slop- 
ing bottom are used to describe the spectral wave components. 
The final chapters conclude with a comparison of the steady 
state solutions produced by the earl ier work and the unsteady 
results obtained here. The numerical results are compared 
with the results for the steady case given by Bowen Cl967j. 
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BACKGROUND 



Changes in mean sea level near a. shorel i ne have been 
studied both theoretically and experimentally. A theoreti- 
cal framework was formulated by Longuet-H i gg i n s and Stewart 
Cl962, 1963, 19643 which dealt with the excess momentum 

flux due to the presence of unsteady wave motion and which 
they termed "radiation stress". Longuet-H i gg i ns and 
Stewart were able to define many of the shoaling effects of 
a train of waves including wave set down and set up utiliz- 
ing this radiation stress concept. This chapter reviews 
the development of the radiation stress tensor and its re- 
lationship to the concepts of set down and set up. It also 
includes a direct approach used by Longuet-H i gg i ns Cl967j to 
derive an expression for set down utilizing the vertical 
momentum equation and the Bernoulli equation. 

A. DEVELOPMENT OF THE RADIATION STRESS TENSOR 

Since the approach of Longuet-H i gg i ns and Stewart is 
rather lengthy and tends to obscure the concepts involved, 
the later and more sysrematic approach by Phi I I ips CI9663 is 
used. The development of the radiation stress tensor and 
the resulting phenomena of set down and set up is kept as 
general as possible. The coordinate system is given by 
Figure (I) where the x-axis is perpend i cu I a r to the shore, 
the y-axis is parallel to the shore, and the z-axis is 
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vertically upward from the still water level. The govern- 
ing equations are the continuity equation and the hori- 
zontal momentum equations. 

The continuity equation is given by 



9p , i , 9pw 
9t 9x . 9z 

I 



0 , 



i = I, 2 



( I ) 



where I, 2 refer to the x, y components. The horizontal 
velocity u. is composed of a mean flow component, U., and 
a fluctuating component representing the deviation from 
mean flow, u.', such that u. = U. + u.'. Since there is 
no mean flow in the vertical, w = w'. Multiplying the con- 
tinuity equation by u. and adding this result to the hori- 
zontal momentum equation 



9 u 

^ W 

produces 



9u . 



9 u . 






> . > 

- + p» ^ 

j 



9 p _ 
^ t R . 
9x . I 

1 



( 2 ) 



9u. 9p u.u. 9p u.w 
I I J I 

0 4 - 

P 9t 9x. 9z 

J 



9x . 

I 



+ R. 

I 



( 3 ) 



Integrating over depth from -h to n , using Leibnitz's rule 
and applying the kinematic free surface and bottom boundary 
condition yields 
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By time averaging this equation term by term and by making 
the following definitions 



M. = J pU.dz = pU. (fi + h) , 

-h 

r 

M . = / pu . ' dz , 

-h 

M. = M. + M. , 

^ M . M . 

U. = 3 —! = U . + ! , 

' p(n + h) ' p(fT+h) 

the following expression is obtained: 



(4) 



9M . 

I 

9t 



+ — { U . M . + 

9x . II 

J 

- J pg(n t h)6.^. } 
= T. t R. . 

I I 




Cpu . 



5 



' u . ' 



J 



i j 



+ p6.^dz 

0 for i r 

1 for i = 



M . M . 

I J 



p ( ri t h ) 



j (5) 

j 



The first term on the left hand side is the local change in 
horizontal momentum flux. In the braces, the first term 
represents the momentum flux produced by the steady state 
flow. The last three terms in the braces contain all the 
unsteady contribution to the momentum flux with "^he hydro- 
static effect subtracted out; this is the momentum flux due 
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to the unsteady motion or the excess momentum flux referred 
to as the radiation stress tensor: 



s,j=y [puj'uj 
-h 



M.M . 

.'+p6..Hdz - ! — ^ 



I J 



p ( n+h ) 



J pg (fi + h ) 6 . j 



( 6 ) 



Equation (5) is then simplified to 



9t 



M. 

I 




{U,M. 



. S,.) 



T. t R. 

I I 



(7) 



T. represents the horizontal force produced by the slope of 
the free surface and is given by 

T. = - pg(n+h) , (8) 



Outside the surf zone it is assumed n << -h and (8) reduces 
to 

T. = - pgh ^ . (9) 

I 



R. is the averaged and integrated frictional stress term. 
These equations are general and apply to all kinds of steady 
and unsteady motion. The only simplifying assumption is 
that mean flow is uniform over depth. The advantage of us- 
ing the "radiation stress" technique to solve physical prob- 
lems is that the second order effects are obtained using 
first order theory. 

B. WAVE SET DOWN 

To present the concept of set down, consider only The x- 
component of horizontal momentum flux equation (7) wi-h 
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waves that propagate shoreward from deep water with their 
crests making an arbitrary angle a with the shoreline. For 
simplicity, the bottom is assumed to be composed of parallel 
contours so that gradients in the y-direction are zero. The 
bottom slope is allowed to vary only gradually so that 
energy reflection from the shoremay be neglected and the 
shoaling effects caused by changes in the bottom can be con- 
sidered in a step-like fashion. By neglecting the frictional 
effects and assuming that any current gradient in the x- 
direction is small, (7) can be written for outside the surf 
zone as 



3M 

X 

3t 




3S 



XX 



9x 



( 10 ) 



Since (10) is an equation in two unknowns (M and n ) / 

X 

a second equation is required for a solution. To provide 
it, the continuity equation, (I), is vertically integrated 
over depth by the use of Leibnitz's rule. This produces the 
conservation of mass flux equation as given by Phillips 
Cl966j where gradients in the y-direction are zero. 



P 




3M 

X 

3x 



0 . 



(II) 



Utilizing I i 
is p roport i ona I 
or to the total 
Long uet-H i gg i n s 



near wave theory, the radiation stress tensor 

2 

to the square of the local wave amp I itude, a , 
energy. Therefore, following the method of 
and Stewart CI9623, the applied force o^ The 
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system travels with the speed of the energy, i.e., at the 
group velocity, c^. Hence the transformation 3/3t + 
c 9/3x = 0 can be applied. Equations (10) and (II) become 
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3M 

X 1 u oO _ 

— * <= 9 ^ 3 ^ ■ 



9S, 

3x 



XX 



( 12 ) 



3M 
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3x 



c Aa 
g 3x 



0 , 



( 13 ) 



for which the solution is 



3x 



or on i nteg ration 



P 



I 



3S 



XX 



, , 2 3x 

(gh - c ) 

^ g 



( 14 ) 



n = p 



I 



(gh - c 2) XX 

g 



( 15 ) 



Since c^ gh in shallow water, equation (15) is a non- 
steady state solution which implies that the mean water 
level increases negatively without bound as the wave moves 
into shallow water. The explanation offered by Longuet- 
Higgins and Stewart Cl962]] for this apparent resonant condi- 
tion is that its effect takes time to build and the energy 
involved is dissipated prior to it reaching any significance 
by the breaking of the wave. 

By imposing steady state conditions and describing the 
unsteady motion using linear wave theory, Longuet-H i gg i ns 
and Stewart Cl9623 subsequently developed the following situ- 
ation for n : 



a^k 



2 s i n h 2 k h 



( 16 ) 
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This Is a second order equation in terms of the local depth, 
amplitude, and wave number. It is apparent that as the wave 
train approaches the point of breaking, q decreases and there 
is a set down of the mean water level. The set down in- 
creases to the point where the wave breaks and other assump- 
tions regarding rj must be appi ied, 

Longuet-H i gg i ns CI9673 also derived this solution in 
another manner without referring to the radiation stress 
term. By integrating the vertical momentum equation over 
depth and time averaging, the total average vertical momentum 
is obtained: 

(p - pw^)^^Q - pgn = 0 . (17) 



A second equation used is Bernoulli's integra 






with the restriction that the flow be i r rota t i ona I . By set- 
ting z = 0 and time averaging, (18) becomes 



I -2-2-2 

p ^ + -^p(u +u +w) ^-)-C = 0, 

'^z=0 2 X y z = 0 ' 



( 19 ) 



where C is at most a constant. From (17) and (19), V 
can be eliminated giving 



n = 



11,-2 

2 g X 



- 2 

u 

y 






C = 0 



( 20 ) 
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From this, the difference in mean sea level is obtained at 
two different points (Xj, 0) and (X 2 > Y 2 ' thereby 

eliminating the constant C, 



An 





) 



z = 0 



] 



1 

2 • 



(21 ) 



The velocities are expressed applying linear wave theory us- 
ing the following relationships: 



^ _ a a ^ ° s a k( z-h ) cos ( kx ' -at+e ) , 

X s I n h kh 

^ _ a o s i n a k ( z-h ) cos ( kx ' -at+s ) ,(22) 

y s I n h kh ' 

w = , u sinh k ( z-h ) s i n ( kx ' -at+£ ) , 

s I n h k h ' 



where x' = x cosa + y s i na and e denotes an arbitrary con- 
stant phase angle. Substituting equations (22) into the 
right hand side of (21) yields 



I 2 ^ - 2 -2, 

( u + u - w ) = 
2g X y 



2 2 
a a 



2 

4g sinh kh 



(23) 



and the difference in mean sea level is then obtained as 



An = - 



a^k 



] ' 

2s i nh 2kh“^2 



(24) 



By assuming the point (X 2 > 'i 2 > infinitely deep 
water so that TI 2 = 0, the value for n at the point (X|, 
0) becomes 



n 



1 a^k 

2 s i n h 2 kh 



(25) 



This solution is very simple and s t ra i gh t f o rwa r d , with the 
same assumptions as before applying to the bottom slope 
and local depth h. This solution is plotted in Figure (2) 
with the strong correlation to laboratory results given by 
Bowen, Innman, and Simmons CI9683. 

As waves progress shoreward they travel from deep water 
to shallow water. The transition in the linear equation 
can be made by utilizing the conservation of energy flux 
(E * Cg = constant). However, when a wave approaches the 
point of breaking, the wave steepens beyond that allowed by 
linear theory and another approximation must be made. It 
is assumed that at breaking 






(26) 



where y is a p report i ona I i ty constant which is generally 
assumed to be equal to 0.78 borrowed from solitary wave 
theory. Then the mean sea level depression at the breaking 
point, rij^, is given by 



Hk = - tV H, 



X 

I 6 



( 27 ) 
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C. SET UP 



Set up is a phenomenon that occurs in the surf zone 
shoreward of the breaker line. Because energy decreases 
shoreward in this zone, a different formulation for the 
wave amplitude must be used than the local value used be- 
fore. It is assumed that the breakers being considered 
are of the spilling type which retain their harmonic charac- 
teristics and gradually decrease in amplitude as they 
progress shoreward such that wave height is defined by 

H = Y(n + h) , (28) 



where y 'S the p report i ona I i ty constant introduced earlier. 
Again by neglecting frictional effects, considering only 
x-direction motion, and by assuming that any current gra- 
dients in this direction are small, equation (7) yields 



3M 3S .- 

X , XX / " , u \ on 

W * 33T- “ ■ ■ 



(29) 



where it cannot be assumed that n << h. In this equation, 
S can be determined by using I inear wave theory and The 

X X 

shallow water approximation giving 



S = 4 E , 
XX 2 ' 



which inside the surf zone is equivalent to 



S 



XX 



_ 3 _ 
I 6 



2 

pgy 



( h + n ) 



2 



By assuming steady state conditions, (29) can be written as 



L- 1 

3x 2 



E 



pg(n + 



h) |I! 

3x 



Since energy decreases shoreward, it is evident that n in- 
creases and produces a set up effect. Solving for n 

n = - kh + C , 



where 




By applying the value of ri|^ at the breaker line, the con- 
stant C can be el iminated and 



n = k(h^-h)+ri[j* (30) 

Where a beach has a constant slope given by h = mx, the set 
up also has a constant slope proportional to the beach slope. 
It should be noted that the simple solution obtained by 
application of Bernoulli's integral outside the surf zone 
is not applicable since the motion is no longer irrotational 



inside the surf zone. 



III. APPLICATION OF A NARROW BAND WAVE 
SPECTRUM TO THE EXISTING THEORIES 



The expressions just considered were made as general as 
possible in order to afford a simple physical understanding. 
Only a monochromatic wave train was considered and only 
steady state solutions were allowed. In this chapter, un- 
steady terms are introduced into these expressions by apply- 
ing I inear wave theory in the form of a simple two frequency 
component wave spectrum simulating a narrow band wave 
spectrum. In this manner an unsteady second order expression 
for the mean water level, q , is derived. Next, using this 
simple spectrum, another solution for n is derived using 
Lon guet-H i gg i ns ' application of the Bernoulli equation. Then 
by i ncorporati ng this spectrum into Iwagaki's first order 
solution for a sloping bottom CI9723 and applying the Longuet 
Higgins' approach, a third solution for q is obtained. 

Final ly, the appi ication of the spectrum to the set up 
phenomenon inside the surf zone is considered. 

A. WAVE SET DOWN 

Consider a simplified wave spectrum, specifically consist 

i ng of two waves which have nearly the same frequency and 

wave number such that 0| ~ ^2 ~ '^l ” ~ 

small. Also for simplicity allow their respective amplitudes 

to be identical and equal to . This amplitude is chosen 

/2 

in order that the variance of the simple wave spectrum and 
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the variance of a monochromatic wave amplitude are the same 
The wave components are to be added linearly sucn that the 
particle velocities become u = U| + U 2 » I 'f the first order 
velocity potential is given locally by 



4 ) 



ag cosh k ( h + z ) 
/2 ck 



s i n ( k . X . 

I I 



at) , i = I, 



where k. is used to indicate direction, i.e. k /k = 

I ' X 

^ A 

and ky/k = sina, then in general u. = - or 



2 

(31 ) 
cosa 



u . 

I 



ag 

/2 c 



^ i cosh k ( h + z ) 
k cosh kh 



CO s ( k . X . - 
I I 



at) . 



(32) 



Since u. = u. + u. in this case, u. becomes 
I I 1 1 2 I 



ag 
/2 < 



i| cosh k|(h+z) 
I cosh k I h 



cos(k. X. - a.t) 
' I ' ' 



k. 



i 2 cos h k 2 ( h + z ) 



cosh k 2 h 



cos ( k . X . - a_t ) . 
I 2 1 2 



(33) 



Since ^1 * *^2 ~ I ~ ^2 ~ subscripts may be 

dropped and (28) can be written as 



ag 



k. 



/2 c 
k . +A k . 






i cosh k(h+z) 



cosh kh 



CO s ( k . X . - at ) 

I I 



+ ( 



I 



I 



k +Ak 



cosh[I( k+Ak) ( h + z)H 
cos h ( k+A k ) h 



cosC( k.+Ak. )x j -(a+Aa) t]} 



(34) 
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This simplifies to 
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/2 c 



k A k 

i cosh k( h+z ) , i Aa 

1 — pr c O S ( — X . - -rs- 

k cosh kh 2 I 2 



t ) cos ( k . X . -at ) 



I I 



( 35 ) 



Similarly, the particle velocity in the z-d i rec t i on , w, is 
obtained and is 



w 



2ag 
/2 c 



s i n h k ( h + z ) 
cosh kh 



Ak. 



cos (■ 



X . 



Aa 

2 



t ) s i n ( k . X . -at ) 

I I 



( 36) 



These expressions state that the particle velocities are 
sinusoidal with a periodically varying amplitude. The modu- 
lating wave is of low frequency and with a wavelength that 
is much longer than those of the original two base waves. 

The unsteady velocities specified by (35) and (36) are 
substituted into (6) and the radiation stress tensor deter- 
mined. The time averaging process is carried out over the 
shorter periods of the base waves. Since the modulating 
wave is of such low frequency, its effects remain in the 
solution. Specifically, equation (6) becomes term by term 



/ 



c k.k. Ak.x. Ak.x. 

^ I r 9 I J ,11 Aa ,, , j j Aa ^ , 

pu.'u.'dz = E — ^ ^ cos( — s — - -p- t)cos( — i ^ t) , 

I J c 2 2 2 2 2 ' 

K 



-h 



( 37 ) 



2 . 



where E = 2 Pga is the total energy density, and c^ and c 
are the group velocity and phase velocity respectively; 



/ 



pdz = pg ( r|th ) + 



o c c « Ak . 

^ ^ - E(^ - 2)cos^(^r^ 

c 



Aa . , 

X|- ^ t) ; 



( 38 ) 



25 







n 




















I 






b‘ 








!r tiffs' r 




and the remaining unspecified term becomes 



M.M . 
I J 



p ( fit h ) 



^2 k. k . 
E ' J 



p ( ri + h ) c‘ 



Ak. 

, I 

cos ( 2 ' ^ j “ 



Ag 

2 



Ak . 

t ) cos (— 2 ^ 



w . - 

J 




(39) 



Substituting (37), (38), and (39) into (6), the radiation 

stress tensor becomes: 



S. . 
I J 



ck.k. Ak. . Ak. 

E — — COS(^7-^ X.- t)C0S(— X . t) t ^ 

c,2 2i2 2j2 2 



Ag 



c Ak . . ^ 

r- / g o \ 2 , I Ag , , 

E(— ^ - 2) COS (^5— X.- t) 

c 2 I 2 



(40) 



p ( n+h ) c 



^2 k. k . Ak. . Ak . 

^ cos(^ X|- — t)C05(-ji .j 2 



-M 

X r - -o- t ; 



Since only second order effects are being considered and 
2 4 

E = 0(a ) the last term may be neglected. It is obvious 
that the proper choice of assumptions will simplify (40) 
immensely, e.g, consider only the x-direction components 
in sha I I ow water, 

k o A k . . _ 

c* 1 AcJ.x ,E 

° — ^ * 2 ■ 

and if the waves approach the shoreline such that the angle 
a = 0, this expression further reduces to 

S = E[l + cos(Ak.x. - Aot) + ^ . (41) 

XX II 2 
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This last expression states that the radiation stress, in 

this circumstance, is composed of a steady stress term and 

/ 

of one that is periodic. 

This derived form of the radiation stress is now used 
to find a solution for ri . Recalling equations (10) and 
(II), and as before neglecting the friction term and assum- 
ing that the current gradients in the x-direction are small, 
the x-component of these expressions become: 



3M 
. > 

3t 



3S, 

3x 



XX 



, 3n 



, 3f! 

— ^ P 3t 



^ = 0 . 



By cross differentiating. 



the M term can be el iminated and 

X 



- P 




+ pg 




afi-, 



2 

3 S 



XX 



3x‘ 



( 42) 



is obtained, 
constant, (42) 




Assuming a flat bottom so that 
becomes 



+ 





XX 



h 



is considered 



(43) 



This is obviously the long wave equation which is being 
forced by the radiation stress term. It has a solution in 
the homogeneous case of 

n = a cos(Ak X - Aat) . 
c X 
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For a particular solution, make the assumption that led to 
equation (41), i.e., consider shallow water and a = 0. Then 

assume a particular solution of the form 



n = A cos(Ak X - Aat)+ B sin(Ak x - Aat) , 

p X X 



By substituting this into equation (43), the solution for A 
and B are obtained 

EAk ^ 



p(Aa^ - ghAk^) 



B = 0 . 



The solution for n is then 



I 2 

j 3 g 

n = Ca - D cos(Ak X - Aat) , (44) 

( - gh) 

Ak^ 

X 

which is a wave with an amplitude that is proportional to 
the steady state solution of Lo ng u et-H i gg i n s and Stewart, 
equation (15). 

If the same assumptions as before are made except tha- 
the bottom is al lowed to vary as h = mx where m is some 
constant slope, a more difficult problem is encountered. 
Equation (42) becomes 



. 2 - . 2 - 9^S 

8r|,^ 9n,„ 9n XX 

-P 4 + pgm 3- + pgmx 4 = ^ — 

9t^ 9x^ 9x^ 



(45) 



This equation again has the characteristics of the long wave 
equation that is being forced by the radiation stress. An 
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analytical solution of this equation would provide a general 
form for n; however, since the depth is al lowed to vary with 
X, the expression for the radiation stress term becomes com- 
plicated and obtaining a particular solution for (45) becomes 
difficult. 

B. WAVE SET DOWN USING THE BERNOULLI INTEGRAL 

Recalling that the second method of Longuet-H i gg i ns to 
determine a solution for n resulted in the equation 

A- I r/- 2 ^ - 2 -2. ^1 

An = T 5 — L(u + u - w ) , 

2g ^ X y z=0 2 ' 



a second unsteady solution may be obtained by applying the 

_ 2 

simple wave spectrum. From equations (35) and (36) u , 

- 2 -2 

u , and w become 

y 

? ? k 9 A k 

- 2 a g / , 2 cosh k(h+z) 2, ^x Aa 

% = — r <— > — XT — <— X - — t) . 

c cosh kh 

? 9 k 9 A k 

- 2 . y^2 cosh k(h+z) 2, Aa 

' c cosh kh 

-2 sinh^k(h + z) 2,^^x ,^*^y Aa j. ^ 

w = g := cos X +— y - ^ t) , 

c cosh^kh z z z 

I o o o _22k„^Ak . 

I /- 2 - 2 -2, _ a g r x,2 2, x Aa 

^ ■ 2g X y ■ z=0 ” 2^2 ^^k 2 x- — 

, , '^y , 2 2/'^y Aa , , 

+ (-r^) COS (— y- -X“ 

k 2 ' 2 



- « A k A k 

- tanh'^kh cos‘^(— ^ x + — 



^ y - t) } . (46) 
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0 and 



Again assuming a = 0 so that k^/k = l> and k^/k 

2 

recalling that c = g/k tanh kh, (46) becomes 




Cl + cosCAk^x - Aat)] . (47) 



n 



2 s i nh 2kh 



This solution is identical to the solution of Longuet-Hig- 
gins (25) except for the unsteady term. The set down 
fluctuates at Aa = Oj - o 2 > which is the difference of the 
radial frequencies of the two component wave spectrum. This 
low frequency oscillation is often described as surf beat. 

The amp I itude of n varies from twice the amp I itude of 
C25) when the two base waves are in phase, to zero when the 
base waves are out of phase. Hence the set down can be con- 
sidered the superposition of a steady component and an 
unsteady component. 

C. WAVE SET DOWN ON A SLOPING BEACH 

In an effort to more accurately represent shoaling wave 
t ra n s f o rma t i o n s , Iwagaki Cl972^ implicitly considers con- 
stant bottom slope resulting in solutions to the long wave 
equation involving the Bessel functions. He obtained first 
order solutions for q and u in terms of the asymptotic form 
of the Bessel and Neumann functions. It seems reasonable 
to assume that by i ncorporat i ng these solutions into the 
Bernoulli integral, a second order solution for q can be ob- 
tained which includes a more realistic rep resen tat i on of the 
bottom effects. Assuming shallow water conditions with the 
depth given by h = mx and progressive waves. 



u 



X 
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t 



n 



a 



( 49) 



/2 



LJq(x'cos 



at + 



N„(X) 



s I n 



at] 



2 . a^x 

where y =4 

gm 

and Neumann functions are: 



The asymptotic expansions of the Besse 



J (u) ~ ^-7- cosCco ~ ^ ~ f 

n f TTW 2 4' 

( 50) 

N ( u ) ~ V““ s i n ( w - ~ • 

n T TTW 2 4 



Substituting these expressions into (48) and simplifying 
yields 



u 



a 

/2 




[s i n (x 



yn 

4 



at ) ] 



(51) 



and from the kinematic free surface boundary condition 



3ri a ^// gmx r • / j.si 

Csin(x - - - ot)] . 



3t 



(52) 



Introducing the wave spectrum as before, squaring and time 
averaging the result, produces 



- 2 

u 

^ z = 0 


a^g/gmx 

2mx^7ra 


Cl 


, ,2Aax 

+ cos ( 

/gmx 


+ Act ) ] 


and 










-2 

w 

Z = f) 


2 2 ^ 

a a /gmx 

2x7Ta 


[1 


, ,2Aax 

+ cost — — 

/gmx 


+ Act)] 



(53) 



( 54 ) 



By assuming that ri 2 is in infinitely deep water as before, 
the second order solution for n is 



D. WAVE SET UP 

Because the nature of the solution for n using linear 
theory was the superposition of a steady state component and 
an unsteady component, it is reasonable to assume that this 
condition persists across the breaker line and the set up 
resulting in the surf zone from the wave spectrum has a 
similar makeup. The steady state component is then derived 
separately in the same manner as in Chapter I I with the 
radiation stress being now defined as the steady state por- 
tion of equation (41), or 



n 




- a^]Cl+cos( 



2Aax 

gmx 



+ Aat)] . (55) 




( 56) 



and aga i n 



n 



s 



kh + C 



( 57) 



where this time the proportionality constant is again 



k 



I t 




Since n. = TZ' YHu. (57) becomes 
D Id d 



n 



s 



k(h. - h) 4- n. . 
D D 



( 53) 
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The unsteady component is assumed to be periodic in 
character; however, unlike the higher frequency waves of 
the train, observations suggest that it is not attenuated 
as it approaches the shorel i ne but is reflected to some 
extent. The formulation is then that of the long wave 
equation as given by Stoker CI9663, 





gh 



8 h 3ri 
3x 3x 



0 , 



( 59) 



and the solution is given by Guza and Bowen Cl9773 as 



n = a{J^(x)sin Act t N^(x)cos Act + 



rCJQ(x)sin(Aat-£)-N^(X)cos(Aat-£)I|} (60) 



where a is the amplitude of the incoming wave (in this case 
a = ri|^), r is the reflection coefficient, and £ a phase 
shift in the reflected wave. The frequency of this wave 
must be the same as that of the modulating wave derived out- 
side the surf zone or Aa. By substituting the asymptotic 
form of the Bessel and Neumann functions, 

{cos ( x~Aat - j) + r cos ( x-Aot-0 ) } , (61) 

TT 

where 9 = +£• Since the set up shoreward was assumed to be 

a composite of the two solutions (58) and (61), 




= k( h -h )+i1l. 
b b 



(I 



■y^Ccos(x-Aat- j) 



■ rcos ( x-Aat-0 ) ] } 



( 62) 
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IV. COMPARISONS 



The results obtained here are compared with the labora- 
tory experiments performed by Bowen, et al.Cl9683; Figure 2. 
They obtained experimental results from a controlled wave 
tank experiment, and in each case only a monochromatic wave 
was considered. Their results compare quite favorably to 
the steady state solutions of Lon g ue t-H i g g i n s and Stewart 
Cl 962]. 

Since it was assumed that the variance of the simple 
wave spectrum was equal to the variance of the monochromatic 
wave, the derived expressions differ from existing theory by 
that amount which is contributed by the modulating wave of 
the unsteady- portion. Comparing Long ue t-H i gg i n s ' solution 
(25) 

- _ _ j_ a^k 
'' 2 s i n h 2kh 



with the unsteady solution (47) 



n 



2 



2 

a 

s i n h 




+ cos ( Ak X - Act ) ] , 

X 



it is found that, by averaging the effect of the unsteady 
component, they are the identical. As mentioned earlier, 
the total effect of the unsteady motion is, in this case, 
to produce a set down that oscillates periodically from 
twice the steady solution to zero times it, i.e., no set 
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down at all. Thus a group of waves which are close in fre- 
quency and wavelength may be said to produce in the near 
shore region a fluctuation in the mean sea level seaward of 
the breaker line. This is also demonstrated by the unsteady 
solution obtained by allowing for a sloping bottom, eq ua t i on 
(.55). It too consists of a steady component and an oscillat- 
ing component. 

The initial solution, equation (44), on the other hand 
is comprised of only an oscillating component. However, this 
solution was derived in a simplistic manner by assuming that 
the depth remained constant and by solving the resulting 
differential equation. While this demonstrates that the 
mean sea level oscillates in a wave-like fashion when driven 
by a group of waves, it does I i tt I e to provide insight into 
the near shore region. 

The set down can be expressed in sha I low water in terms 
of deep water conditions by using conservation of energy 
flux. 



E c 




con stant 



where 
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The subscript o refers to deep water and the subscript s 
refers to shal low water. The set down in shal low water 
by equation (25) is then given by 



n 



2 1/2. ,-3/2 

a g (mx) 
o^ 

40 



( 63) 



where in shallow water sinh 2kh ^ 2kh and h = mx for a con- 
stant sloping bottom. 

The set down derived from I inear theory can be compared 
to that derived from linear theory for a sloping bottom, 
equation (55), by also expressing the shallow water amplitude 
for sloping bottom solutions in terms of deep water condi- 
tions. This is accomplished by patching Stokes solution 
at the off shore to the shallow ater solution. For a smooth 
match it is required that, from Friedrich [19481, 

2 

2 _ ^o ^ 

^ s 2m * 



The steady state solution of (55) can now be expressed in 
terms of deep water conditions 



n = - 



4a 



2 . 2 

o 1/2. ,-3/2 ^ ° , ,-1/2 

g ( mx ) + (gmx) 



2 

9 0 I / 7 

= n I + ( gmx ) " 



(64) 



Hence, the sloping bottom solution results in slightly 
greater set down as compared with I inear theory. Numerical 
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figure 2' Profile of Mean Sea Level 



results of these equations, (63) and (64), are plotted in 
Figure 2 . 

Close to the break point, most theoretical results fai 1 
to compare favorably with the experimental results. Theory 
increases set down rapidly near the break point while ex- 
perimental results tend to flatten out. This is due most 
likely to the failure of linear theory to adequately describe 
b rea k i ng waves . 

Although simplistic in form, the application of the two 
frequency wave spectrum demonstrated that fluctuating 
values are obtained for the set down and set up phenomena, 
i.e., time dependent solutions were obtained for these con- 
cepts. Hence, a group of waves, similar in frequency can 
be expected to produce a periodic variation in the mean sea 
level which becomes most apparent when they encounter a 
beach . 
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V. CONCLUSIONS 



In this study, unsteady solutions for set down and set 
up were derived using a simple two component wave spectrum. 
The spectral wave components were described using I inear 
wave theory for both horizontal and sloping bottoms. The 
sloping bottom solution involving the Bessel functions gave 
a slightly greater set down compared with that of the 
linear wave theory solution. However, both unsteady solu- 
tions for set down showed, that to at least a first approxi- 
mation, a steady and a fluctuating component, the steady 
component being identical to the earlier steady state 
so I u t i on s . 

The set down at the breaker line acts as the boundary 
condition driving the set up inside the surf zone. The set 
down solution showed that the steady and unsteady components 
are simply additive. Hence, inside the surf zone a solu- 
tion is composed of the steady set up and a long wave, i.e., 
surf beat, which is driven by the fluctuating condition at 
the breaker line. 
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